Core-crust transition pressure evolution in post-glitch epoch for a Vela-type pulsar 
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We propose that the post-glitch epoch in a Vela-type pulsar corresponds to a transition time in 
which the broken crust is readjusting its temperature and core- crust transition pressure after the 
changes produced by the glitch. To describe the evolution of the pulsar in this epoch we use a 
sequence of stationary and axisymmetric relativistic models of slowly rotating neutron stars with a 
surface layer crust. Previously, we have formulated the matching conditions on the surface of the 
star in terms of physical properties (total mass of the star, core-crust transition pressure,...). 

PACS numbers: 97.60.Gb, 97.60.Jd, 04.40.Dg, 95.30.Sf 



Many pulsars show sudden spin jumps, glitches, which 
are superimposed to the gradual spin down due to the 
continued loss of angular momentum suffered by the star 
[H . After a relaxation period of the order of hundreds of 
days (post-glitch epoch), the star almost recovers its usual 
rhythm of spin down. The study of the properties of both 
glitch and post-glitch epoch are essential to understand 
the structure of the neutron star that models the pulsar 
1- 

Glitches have been observed in most of the younger 
isolated pulsars. Particularly, more than a dozen glitches 
with period variation of the order of 1CP 8 and 1CP 6 have 
been observed in the Crab and Vela pulsars respectively 
since their discovery. Many mechanisms have beenpro- 
posed to explain glitches [3[ and post-glitch epoch 

Most of the models that deal with glitches are asso- 
ciated with the layer structure that neutron stars are 
thought to posses. Basically two regions of the star in- 
terior can be differentiated: the core and the crust. The 
properties of the crust are very different to those of the 
core, and it is thought that this region has a solid crys- 
talline structure similar to a metal [5(. 

Let us mention that in the thermal glitch model pro- 
posed by Link and Epstein the glitch is explained by 
a sudden deposition of thermal energy in the inner crust, 
energy that would come from a rupture of the solid crust. 
The inner crust superfluid is very sensitive to star quake 
generated temperature perturbations, affecting to the ro- 
tation behavior of the star. Furthermore, it is notewor- 
thy that the dynamical properties of the neutron star 
depend strongly on the transition pressure between the 
star's core and the crust as have been pointed up by 
Cheng, Yuan, and Zhang ;d] and Lattimer and Prakash 
0. 

According with those two ideas, we propose that the 
post-glitch epoch can be understood as a transition time 
in which the broken crust is readjusting its temperature. 
The temperature excess that the inner layers of the crust 
has just after a glitch, would gradually dissipate during 
the post-glitch epoch until the crust reaches its equilib- 
rium temperature. Our hypothesis is that the cooling of 



the exterior layers of the star would lead to a sensitive 
variation of the transition pressure between the core and 
the crust. In consequence, the post-glitch epoch will cor- 
respond to a readjustment epoch for this transition pres- 
sure. Because many dynamical properties of the neutron 
star strongly depend on the core- crust transition pres- 
sure, the gradual cooling down of the inner layers of the 
crust would explain its evolution all along the post- glitch 
epoch. 

In order to be able to control the core- crust transition 
pressure, we will consider that the crust is thin enough 
so that we can treat it as a surface layer that envelops 
the star's core. We will consider the case of a relativistic 
star with rigid and slow rotation, so we can make use of 
Hartle and Thorne formalism [§]. The introduction of a 
surface layer of energy density will change the boundary 
conditions on the surface of the star. We will show that 
it is interesting to use the total mass, the central density 
and the core- crust transition pressure as the parameters 
that characterize a configuration. 

When performing the simulation of the post-glitch 
epoch, we must take into account that there are two very 
different time scales in the evolution of the star: the post- 
glitch epoch duration, that is of the order of hundreds 
of days, in contrast with the rotation period of the star 
that is of the order of seconds or even milliseconds. This 
important feature allows us to approximate each evolu- 
tion state of the star during the post-glitch epoch, using 
an hydrodynamic equilibrium configuration with the ap- 
propriate parameters (total mass, central density, and 
core- crust transition pressure). We will fit the evolution 
of the angular velocity of the star with a sequence of 
configurations of different cove-crust transition pressure, 
assuming that neither the mass nor the central density 
changes during the post-glitch epoch. 

Model of a slowly rotating neutron star with a 
surface layer crust 

We assume that the neutron star is in permanent rigid 
rotation. The star rotates with constant angular veloc- 
ity Q around an axis, so the resulting space-time is ax- 
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isymmetric and stationary. Appropriate coordinates can 
be chosen so that the metric can be written as follows 
ds 2 = -H 2 dt 2 + Q 2 dr 2 + r 2 K 2 [d9 2 - sin 2 9(d<f> - Ldt) 2 ] . 
We assume that the rotation of the star is sufficiently slow 
and, hence, the Hartle-Thorne perturbative solution can 
be used Q. To second order in the angular velocity of 
the star f2, it can be proved that the components of the 
metric have the following expression: 
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We have to solve the Einstein equations in the interior of 
the star, where the matter is described by a perfect fluid 
tensor T^ v = (p + p)u l± u v + pg^v with equation of state 
p = p(p) (here is the four-velocity, p is the proper 
energy density, and p is the pressure of the fluid) and in 
the outer region, where space is empty. Also, we have to 
impose appropriate boundary conditions in the surface 
of the star. 

It is very desirable to introduce a new radial coor- 
dinate R. We can define a coordinate transformation 
r -> r(R,9) = R + £(R,0) 0, so that surfaces of 
R = const, has constant density, i.e., this coordinate is 
adapted to the surface of the star. 

We assume that the core of the star is surrounded by a 
thin surface layer (simplified crust model). Then we find 
that in the surface of the star R = a, the stress-energy 
tensor is non null and it can be written as follows: 



Tr(R) = Pc(9)u^ c (a)6(R-a) 



(2) 



Hence, we also assume that the matter of the surface layer 
can be described as a perfect fluid. We allow this fluid 
to rotate with its own angular velocity f2 c . For the sur- 
face energy densities we will consider, we can despise the 
contribution of the pressure to the surface stress-energy 
tensor. Therefore the surface energy density can be writ- 
ten to second order as follows: 



p c (6) =e + Se(9) = e + Se + 6e 2 P 2 (9) 



(3) 



where e, Seo, and 8e 2 are constants. 

We must impose certain junction conditions on the in- 
terior and exterior solutions taking into account the in- 
troduction of a surface layer of matter on the border of 
the star. We use and intrinsic formulation of these con- 
ditions: the first fundamental form of the hypersurface 
R = a has to be continuous and the second fundamental 
form has to have a jump given by the stress-energy tensor 
on the surface of the star Q . Expanding these junction 
conditions in terms of f2, we find the conditions for each 
one of the metric functions 



- Distribution of pressure inside the star 
p (R) for different transition pressures 



FIG. 1: The intersection between Pi nt , for different transition 
pressures, and p(R) determines the radius of the star 



Order zero: we obtain important conditions for the 
mass function and for the pressure of transition between 
core and crust 
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The equation (|4a[) can be interpreted as giving the total 
mass of the star to zero order (M ext ), as the addition of 
three terms, the first one representing the interior mass 
(Mi n t = 47r J® pR 2 dR), the second one representing the 
mass of the crust and the third one being a negative 
bounding energy term. 

The equation (I4bjl shows that the introduction of a sur- 
face layer of energy density causes a discontinuity of the 
pressure at the border of the star. In the inner part of the 
surface the pressure is pi„t and outside the star the pres- 
sure is null. In our model, pi n t represents the core- crust 
transition pressure. We interpret (|4b|) as an equation 
giving the core- crust transition pressure in terms of the 
total mass of the star, the interior mass, and e. 

This condition is used to obtain the radius of the star, 
as can be seen in fig. [TJ from (|4b[) we define several 
functions Pint(R), for different transition pressures, and 
we study the intersection between these functions and 
the distribution of pressures in the interior of the star. 

First order: we obtain continuity of the rate of 
rotation of the inertial frames, and an expression for the 
total angular momentum: 

u{a) ext = w(a)i n t = w(a) (5a) 

J = ye A < a '»' (e-M«)*»« [d R u(a)] mt + 16neU^5b) 

where ~ c = f2 c — uj(a). 

It can be seen [8( that the first term of expression 
(f5b|) has the same sign as f2 (in our case always pos- 
itive). However, the second term depends on the sign 
of £l c . Hence, contra-rotating configuration are possi- 
ble (the core and the crust contra-rotate), and even, we 
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could find certain critical configuration in which the total 
angular momentum is null. 

Second order: we obtain the following conditions. 

A[r*M)] = (6) 

where r* = r + £ + (£2 + r k 2 )P2(0) is defined as in 0, 
so both the mean radius f* — r + £0 and the eccentricity 
e = [—3(^2 + ^/r)] 1 ! 2 are continuous. We also obtain 
the next conditions 

A [e A ^m (a)] = 47ra 2 fe + ^ a 4 e -^( Q ) £ ^ 2 + 

- i£ (a)A[3 e - A W+ e A ( a )(8™ 2 p(a)-l)] (7a) 
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+ ieo(a)A[F 3 (a)] = -^a 4 e^) £ ^ 2 (7d) 

where we have defined A[/(a)] = /|s ea!( — /|n <n , and 

Fi(a) = e 2A(a) (l + 87ra 2 p(a)) - 1 (8a) 

F 2 (o) = e A(a) (7 + 87ra 2 p(a)) -3e~ A ( Q ) (8b) 

F 3 (a) = 64 7 r 2 a 4 e 3A(a V(a)p(a) + 87ra 2 e 3A(a) ( / o(a) -p(o)) 

+ 87ra 2 e A ( a )(3p(a)-p(a))-e 3A(a) -3e- A ( a ) (8c) 

Basically, this conditions fix the jump in the second order 
perturbation of the mass function, and also fix the second 
order perturbation on the surface density (5eq and & 2 ). 
Finally the last of the second order conditions fixes the 
value of the angular velocity of the crust: 

n c = w(a) ± (O - u(a)) (9) 

We obtain two possible configurations for the same star 
core and transition pressure: one with the crust co- 
rotating with the interior fluid and with identical ve- 
locity, and a special configuration with contra-rotating 
crust. Current research on this strange contra-rotating 
configurations is under development, and it will be pre- 
sented elsewhere. In what follows in this paper we will 
use the co-rotating configuration. 

In a summary, to obtain a model of a neutron star with 
a surface layer crust in permanent slow rotation, for a 
given equation of state of the core, we have to fix three 
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FIG. 2: Mass and radius of realistic neutron stars. 

parameters. We choose the total mass of the star, the 
central density, and the core- crust transition pressure. 
Then, we integrate the Einstein equations inside the star 
(using a realistic equation of state for the matter) and 
we match this solution with the exterior solution found 
in Q using the junction conditions explained above. In 
fig. [2l we show multiple configurations for both static 
stars and rotating stars, with a crust of approximately 
10% of the total mass. The angular velocity displayed is 
f2| = M ^ t (mass shedding angular velocity [8|). On the 
static configurations without crust we write the value 
logl0(p c ). We have used an equation of state for high 
pressures due to Glendenning [n| (K=240 MeV, s£ = 
0.78 and x a = 0.6) and the BPS equation of state for the 
lower pressures. Similar results have been obtained for 
the analitical fit of the Sly equation of state presented in 
1- 

Post-glitch epoch of a Vela-type pulsar 

Now we are going to study the post-glitch epoch of 
a Vela-type pulsar. The typical evolution of the angu- 
lar velocity for a Vela-type pulsar during__gZi£c/i time and 
post- glitch epoch is presented in fig. [3] 11]. It is interest- 
ing to note that in this problem there are two different 
scales of time, one related with the rotation period of the 
pulsar, and the other one with the duration of the post- 
glitch epoch. This fact allow us to describe the evolution 
of the pulsar in the post-glitch epoch using a sequence 
of stationary and axisymmetric configurations with the 
appropriate parameters to follow the evolution of the an- 
gular velocity of the pulsar. We assume that during the 
post-glitch epoch neither the total mass nor, the central 
density of the star are changed. Hence, in our descrip- 
tion, the changes in the pulsar during this epoch are due 
to readjustment of the core- crust transition pressure. 

In order to control the core- crust transition pressure, 
we assume that the crust is thin enough to be treated 
as a surface layer crust. Then to obtain the sequence 
of stationary and axisymmetric configurations mentioned 
above, we use our model of slowly rotating neutron star 
with superficial crust. 



4 




-100 -50 50 100 150 200 



t(d) 

FIG. 3: Evolution of on a typical Vela pulsar glitch. The 
points correspond to the sequence of models. 
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FIG. 4: Evolution of the core- crust transition pressure and 
eccentricity during the post-glitch epoch. 

Let us describe our analysis. We assume that the 
neutron star has 1,44m© and central density 1.279 ■ 
10 15 gcm -3 . At the beginning of the post-glitch epoch this 
star has a crust with the 10% of the total mass. Since 
then, we adjust the transition pressure to fit the angular 
velocity of the model to the evolution shown in fig. [3] 
The points in this figure correspond to the elements of 
the sequence of models used to describe the evolution of 
the pulsar in the post-glitch epoch. 

We present our results in fig. [4] that shows the evo- 
lution of the core- crust transition pressure and eccen- 
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FIG. 5: Evolution of the surface matter distribution and 
quadrupolar moment during the post-glitch epoch. 



tricity, and fig. [5] that shows the evolution of the sur- 
face matter distribution and quadrupolar moment. We 
obtain that the core- crust transition pressure decreases 
with time and causes a progressive reduction in the ec- 
centricity, and a movement of matter on the surface of the 
star, from the equator to the poles. Also, the quadrupo- 
lar moment decreases in time all along the post-glitch. 

In consequence, we think that the post-glitch epoch 
could be understood as a readjustment epoch for the 
core- crust transition pressure. Our conjecture is that 
the temperature excess that the inner layers of the crust 
have just after a glitch, would gradually dissipate during 
post-glitch epoch until the crust reaches its equilibrium 
temperature. These changes of temperature produce the 
evolution of the core- crust transition pressure and the 
rest of the properties of the star. 
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